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Abstract 

Charged pion spectra measured in 58 Ni- 58 Ni collisions at 1.06, 1.45 and 1.93 
AGeV are interpreted in terms of a thermal model including the decay of A 
resonances. The transverse momentum spectra of pions are well reproduced 
by adding the pions originating from the A -resonance decay to the component 
of thermal pions, deduced from the high transverse momentum part of the 
pion spectra. About 10 and 18% of the nucleons are excited to A states at 
freeze-out for beam energies of 1 and 2 AGeV, respectively. 

PACS numbers: 25.75.-q, 25.75.Dw 
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Relativistic heavy ion collisions are a unique tool to investigate hot and dense nuclear 
matter in the laboratory In these studies the measurement of pions has attracted 

much interest since they can be regarded as a potential probe of the processes beyond mere 
nucleon-nucleon binary interactions in heavy ion collisions • According to the results of 
microscopic model calculations pions are, in principle, expected to contain information 
from the early stage of the fireball to the freeze-out, since they are produced continuously 
during the whole reaction process. However, the large pion-nucleon cross section destroys the 
memory to early times, and the signatures of the reaction dynamics have disappeared by the 
time of detection. Therefore one will, most likely, observe only the conditions at freeze-out. 
On the other hand, this knowledge is important in itself, because it gives information on the 
degree of equilibration. To a large extent, a model that assumes thermal and hadrochemical 
equilibrium can explain the various particle yield ratios including multistrange particles at 
AGS and SPS energies @. 

The possibility of exciting A(1232) resonances (called A in the following) and hence to 
form so-called A matter at incident energies between 1 and 10 AGeV has been investigated 
theoretically in RQMD simulations |J. The model predicts about 10 % of the participant 
nucleons being excited to A during the hot and dense phase at 1 AGeV, while at 10 AGeV 
more than 50 % of the nucleons are excited to baryonic resonances, of which only about half 
are present as A . Pions from the resonance decays show a different phase space distribution 



depending on their decay momentum [10|, and especially the A - decay kinamatics has been 



made responsible for the low transverse momentum (p t ) enhancement at incident beam en- 



ergies from 1 to 15 AGeV [nj-1^]. The concave shape found in the Boltzmann (or invariant) 
spectra of pions at beam energies from 0.5 to 200 AGeV []T~3| pT8f| confirms the importance of 
the resonance decay. In this paper we will investigate the A population at freeze-out, using 
the spectral shape of pion distributions obtained in the Ni+Ni system at incident energies 
varying from 1 to 2 AGeV. 

The experiment has been performed with the FOPI detection system [ fr9| , |2"0" | at the 
heavy ion synchrotron SIS at GSI by bombarding a 58 Ni target of 225 mg/cm 2 thickness 



(~ 1 % nuclear interaction length) with 58 Ni beams of 1.06, 1.45 and 1.93 AGeV. For the 
analysis presented here the central drift chamber (CDC), working in the 0.6 T field of a 
solenoidal magnet has been used for particle identification, and the forward plastic wall for 
the centrality determination. This azimuthally symmetric wall covers polar angles Oiab from 
7 to 30 degrees, measuring energy loss dE/dx and time of flight of the charged fragments 
the multiplicity of which (PMUL) is used for centrality selection. The CDC covers 9 Lab from 
30 to 150 degrees; identification of the pions and other light charged fragments is achieved 
via dE/dx and momentum measurement. To compare the rapidity distributions at different 
beam energies we define the normalized rapidity y^ = y/y cm — 1, where y cm is the rapidity 
of the center of mass system (cm.). Throughout this paper we use the convention % = c = 
1. 

We have previously observed that the Boltzmann spectra, d 2 N/m 2 dm t dy^ with m t = 



p\ + m 2 being the transverse mass and the pion mass, of ir^ measured in reactions 
between 1 and 2 AGeV show a concave shape which can be fitted by the sum of two expo- 
nential functions with a x 2 P er degree of freedom ~ 1 in all rapidity bins |T5 |. Fig. [I] shows, 



as example, the Boltzmann spectra of tt~ at 1.93 AGeV with a cut in PMUL selecting the 
uppermost 100 mb (or 4% of the total reaction cross section, called PM100 in the following) 
for different rapidity windows. Only the statistical errors are given in the spectra; the sys- 



tematic error in the absolute normalization [fTq| , estimated to about 10 %, is not shown. It 



has been demonstrated that at freeze-out the A -decay pions contribute dominantly to the 
low kinetic energy part, whereas thermal pions contribute predominantly to the high kinetic 
energy part of the spectra both at SIS and AGS energies |TT]-0. Therefore we attribute 
the high-pt component of the experimental spectra as being due to thermal pions. 

In order to decompose the pion spectra quantitatively we calculate the part arising from 
A decay by a Monte Carlo routine. For the phase space distribution of the A at freeze-out, 



we use the isotropically expanding blast model proposed by Siemens and Rasmussen |21 
d 2 N A r°° , „ . , „ ^ Elffr , T.sinha T 



r°° I sin h ry 7 

, , . dm A P A (m A )E A e-^ T [( lr + —) — cosh a], (1) 

p t dp t dy Jm N +m-K E A a E A 



where is the nucleon mass, E A = Jp\ + m A andpA are the total energy and momentum 



of the A in the cm., 7 r = 1/y 1 — (3%, and a = ( r y r /3 r pA)/T. For the freeze-out temperature 
(T) and radial flow velocity (/3 r ), the results from Ref. [|15j are used, e.g. T = 92 MeV 
and /3 r = 0.32 at 1.93 AGeV, obtained in analyses of the high-p t part of pion, proton and 
deuteron spectra at midrapidity. We adopt an isotropically emitting source scenario that 
describes the high p t part of the pion spectrum very well. It should be noted, however, that 
the proton rapidity distributions are more elongated and the apparent temperatures near 
target rapidity are smaller as compared to the expectations from an isotropically radiating 
|15|| . Most probably this is caused by the interaction of participating nucleons with 



source 



some spectator remnants or partial transparency of the system. For those scenarios a reliable 
prescription for the A phase space distribution does not exist at the moment. The A mass 



distribution, -Pa(^a), is parameterized by a relativistic Breit-Wigner form [22] as follows: 



1 T 2 (q) 

PA(m A ) oc (_)__^_^_, (2) 

where 

q 1 + (Riq) 2 + {R2qr 
is the momentum dependent width of the A , q is the momentum of the pion in the A rest 
frame. We chose f A = 112 MeV, q = 227 MeV, i?i = 0.0042 MeV" 1 and R 2 = 0.0032 
MeV -1 . The maximum of P A {m A ) is fixed to the free A mass, i.e. m A = 1232 MeV. 

After superposition with the spectra of thermal pions the results of the model calcu- 
lation are shown by the black solid lines in Fig. [I]. We use for all rapidity bins only one 
normalization factor for the relative abundance of the decay pions with respect to thermal 
pions. The overall agreement between the model and data is rather good except at very low 
m t (m t — m n < 0.12 GeV), where the data exceed the model calculations at all rapidities. 
In order to explain this difference, one may need to consider the decrease of the effective A 
mass as shown, e.g. by the thermal model including the 7r— nucleon loop correction to the 



A self-energy calculation ||12|| . The decomposition of it spectra into ttt (dashed line) and 



7Ta (dotted line) at midrapidity is also displayed in Fig. [I], demonstrating that the \ow-p t 
enhancement of the pion spectra is consistent with the contribution from A decay The 
influence of the p t resolution (<J Pt /pt) on the 7Ta spectra which worsens from 4 % at p t < 0.5 
GeV to about 9 % at p t = 1.5 GeV |23| has been tested, too: a slight effect can be found 



only for m t — > 1 GeV, and it does not affect our results. 

The influence of higher-mass baryonic resonances (e.g. N*(1440), N*(1520) and N*(1535)) 
has been studied using the hadrochemical equilibrium model [§,H]. The model parameters, 
chemical freeze-out temperature Tc and baryon chemical potential are adjusted to re- 
produce the measured particle yields of high-p 4 (thermal) pions, protons, deuterons and rj's 
|T5|,|25[|. The temperatures Tq from the hadrochemical equilibrium model agree with those 



from the radial flow analysis, T, within 8 % p4 |. The Boltzmann spectra of pions from 
the N* resonance decay at midrapidity are also shown in Fig. p], calculated with the proper 
branching ratios under the assumption of isotropic expanding thermal distributions of pri- 
mary resonances (Eq. [I]) with a fixed mass. The total contribution of N* resonances to our 
estimation of the A abundance is at most 5%, in agreement with earlier RQMD and BUU 



studies PJ2^|; hence this contribution is neglected in the current analysis. 



The concave shape in the pion spectra can also be found in transport models [2B,27|, 
although the high-p t pions are of different origin there. Within these models, all pions stem 
from the decay of resonances. Pions in the high-pt tail have decoupled earlier from the system 
(free pions), when more energetic resonances were populated. The pion spectra are obtained 
by integrating over the complete compression expansion cycle. However, when analysing 
the result of the reaction under the assumption of a sudden freeze-out, the spectrum of the 
free pions exhibits an exponential shape as obtained for the thermal pions in the present 
equilibrium model. 

To estimate the yields of all pions (jiau) an d that of the \ow-p t pions (j[J jOW ) we have 
integrated the fitted function (sum of two exponentials) over the whole range pt = to oo 
and for the \ow-p t exponential only, respectively. Since the slope parameter of the low p t 
component from the two exponential fit agree with those of the pion spectra from A decay 



within 6% for all beam energies, we simply take tti ow multiplied with the proper isobar ratio 
as measure of the total number of A in the system. The resulting dN/dy^ spectra are 
shown in the top panel of Fig. ^| for the three beam energies under the PM100 cut. Note 
that only data for y^ < are measured, and then reflected to y^ > 0, employing the 
symmetry of the colliding system. The total 7r~ yield per event, n(ir~), is summarized in 
the second column of Table |. The number of participant nucleons (A part ) is calculated via 



the geometrical model ||23|| ; for PM100 events in Ni + Ni collisions one obtains 101 ± 8 
nucleons in the central fireball. In comparison to the earlier BEVALAC data for Ar + KC1 
@, our measurement of the 7r~ multiplicity per participant nucleon, n(7r~)/A part , is about 
15 % lower than the mean value given there; within the error bars both measurements are 
compatible. The ratio ^Zow/^au i s shown in the bottom panel of Fig. |2|, where the dashed 
lines represent the mean of the ratio for the whole rapidity bins given also in the third 
column of table I. The systematic error of ^low/^au introduced by using the exponential fit 
at low p t rather than the 7Ta spectrum is estimated to less than 2%. The ratio decreases 
from about 77 % to 66 % as we increase the beam energy from 1.06 to 1.93 AGeV. At 14.6 
AGeV this ttlowI^au ratio was determined to 33 ± 5 % JT5| |. 

Having this information, we can now estimate the A abundance n(A) at freeze-out by 
means of the following relation: 

n(A) = n(7r-) x f tsobar x (%^), (4) 

where 

= n(7T+ + 7T° + 7T-) _ 6 (Z 2 + N 2 + NZ) 
isobar- - 5 N 2 + NZ ^> 

is the prediction of the isobar model |1]] (N is the number of neutrons, Z the number of 
protons). For the energy range considered here, this model describes the relative abundance 



of pions in central collisions within 10% [P3,Z^|. For Ni one obtains fi SO b ar =2.8A. Note 
that for the present system iV is close to Z, hence one does not expect big deviations from 
fisobar=3 in any model. Eq. f| contains the implicit assumption that the ratio ttlow/^au is 



the same for tt + , tt° and tt~ . From our analysis the spectra of 7r + and 7r~ at a given rapidity 
are known to be identical for transverse momenta from 170 MeV to about 800 MeV, the 
value up to which 7r + can reliably be identified in the experiment. Below pt = 170 MeV, 
the 7r~ yield starts to exceed that of the 7r + , which is attributed to the Coulomb effect. The 
results for the ratio n(A)/n(nucleon + A) are summarized in the last column of Table |. 
The number of nucleons excited to the A resonance at freeze-out increases from 10 ± 2 % 
at 1.06 AGeV to 18 ± 4 % at 1.93 AGeV. These values are diplayed in Fig. || together with 
the result of 36 ± 5 % from a measurement at 14.6 AGeV []13| |. Having determined the ratio 
n(A)/n(nucleon + A), one can estimate the freeze-out temperature within the context of 
the hadrochemical equilibrium model || . The results are summarized in Table [II] together 
with those of the radial flow analysis jUj; both analyses are in good agreement. 

In summary, we have studied the pion phase space distribution in terms of the thermal 
model including pions from the decay of A resonances in addition to the thermal pions. The 
good agreement between the data and model gives confidence that the \ow-p t enhancement 
in the pion spectra is due to the pions from the A decay. When the beam energy is increased 
from 1.06 to 1.93 AGeV this \ow-p t enhanced part of the pion spectra decreases from 77 % to 
66 % of the total pion multiplicity, while the number of nucleons excited to A at freeze-out 
increases from 10 to 18 % of all nucleons. The freeze-out temperature which is estimated 
from the ratio n(A)/n(nucleon + A) is found to be in good agreement with the results of 
the radial flow analyses. The direct reconstruction of the A by an invariant mass analysis 
of (np)- pairs is underway which will allow to compare the results of the two different 
analyses in the future. 
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TABLES 

TABLE I. Summary of results for Ni + Ni collisions at the various incident energies (centrality 
cut PM100): Average number of tt~ per event (2nd column), fraction of the low- 71^ pion component 
(3rd column) and derived relative number of nucleons excited to A at freeze-out (last column). 

E beam /A(GeV) n(TT-) *I ow /*m(%) n(A)/n(nucleon + A)(%) 

1.06 3.6 ± 0.4 77 ± 4 10 ± 2 

1.45 5.8 ± 0.6 68 ± 4 13 ± 3 

1.93 8.5 ± 0.9 66 ± 3 18 ± 4 



TABLE II. Comparison of T in MeV obtained in this work (A) and in a radial flow 
analysis |I3 . 



E beam /A(GeV) A Flow 

1.06 75 ± 5 79 ± 10 

1.45 80 ± 7 84 ± 10 

1.93 89 ± 9 92 ± 12 
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FIGURES 

FIG. 1. Boltzmann spectra of 7r~ for different windows in normalized rapidity (Ni + Ni 
collisions at 1.93 AGeV under selectivity cut PM100, see text for definition); the two upper curves 
have been multiplied by the factors 10 2 and 10 4 , respectively. The black solid lines represent the 
model calculations including the thermal pions (dashed) as well as pions from the A decay (dotted) . 
Pions from the decay of higher baryonic resonances, e.g. from iV*(1440) (dash-dotted), iV*(1520) 
(grey solid-a) and iV*(1535) (grey solid-b) decays, are shown as well. 



Fig. 1 
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FIG. 2. Top: Rapidity spectra of ir~ at 1.06, 1.45 and 1.93 AGeV under the centrality 
cut PM100. The solid circles represent the integrations from p t = to oo of the fit with two 
exponentials (tt^ u ), while the open circles represent the integrations of the low-pj exponential fit 
only (iTi ow )- Bottom: Ratio of ^Zow/^AU' * ne dashed lines represent the mean values of the ratio 
for all rapidities. 
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FIG. 3. Ratio n(A)/n(nucleon + A) at freeze-out as a function of the beam energy. 
The open square is the result of the analysis of experiment E814 at 14.6 AGeV [ITS 



Fig. 3 
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